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Abstract 

We consider the implications of the recent determination of the universal in- 
frared background for the propagation of photons up to 20 TeV from the active 
galaxy Markarian 501 as observed by HEGRA. At 20 TeV the mean free path for 
photon-photon collisions on the infrared background would be much shorter than 
the distance to Markarian 501, implying absorption factors of the order of exp(— 10), 
or greater, and consequently an excessive power output for this active galaxy. Pos- 
sible solutions of this problem are discussed. 
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1 Introduction 



The general form of the spectrum of the universal infrared background radiation was 
predicted by Primack et al. [|J, Malkan and Stecker 0, Dwek et al. and Fall et 
al. [[§]. The origin is currently thought to be due to star formation in the early universe 
producing starlight predominantly in the 1 micron range which was then re-processed into 
the 100 micron range due to heating of surrounding dust. Early data from the DIRBE 
and FIRAS instruments on COBE (see ref. for a review) and from ISO || were mainly 
upper or lower limits which did not severely constrain the models. New results from these 
experiments are becoming available due to a more sophisticated treatment of the galactic 
and solar-system radiation foreground subtraction. In particular, new determinations of 
the infrared background intensity have very recently become available at 1.25, 2.2 and 
3.5 microns |J, 15.0 microns J7|, 60 and 100 microns ||, and 140 and 240 microns §J. In 
Fig. 1 we show the predictions of Malkan & Stecker || together with the new observations 
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and recent lower limits |TIJ based on infrared galaxy source counts and a model including 
the effects of an infrared burst phase of ultraluminous infrared galaxies. The observations 
are in excellent agreement with prediction inasmuch as the expected trend (two peaks 
and valley) is found, the agreement being particularly good between 1.25 and 15 microns. 
However, the level of the observed flux in the 60-240 micron range is higher, and in the 
case that this is a truly universal background it has important consequences both for 
star formation in the early universe, and for intergalactic absorption of TeV gamma-rays 
which leads to a crisis in our understanding of high energy phenomena. 

Immediately following the discovery of Markarian 421 as a TeV source |l]J it became 
clear that interactions with a universal infrared background may have a strong influence on 
the propagation of TeV gamma-rays through intergalactic space |12 , |13| , in the same way 



that the microwave background affects the propagation of 1000 TeV gamma-rays fll~4| , |15j . 
The observation by several telescopes of a very high level of emission by Markarian 501 
during 1997 (see ref. |TB| for a review) showed that the spectrum extended to 10 TeV 

that if 
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and beyond. At that time, it was pointed out by Stanev and Franceschini 
the spectrum continued much beyond 10 TeV that we would "have to revise our concepts 
about the propagation of TeV gamma-rays in the intergalactic space, and that something 
complicates the process." This has now come about as the analysis of HEGRA data for 
the whole of 1997 shows that the spectrum of Markarian 501 extends well beyond 10 TeV 

m. 



2 Implications of the new IR observations for the 
propagation of Gamma-rays through the universe 

Photon-photon pair production is a very elementary process in QED, and can safely be 
calculated at energies relevant to TeV astronomy. For propagation of energetic parti- 
cles (gamma-rays) of energy E y , mass m 7 (=0) and velocity /3 7 c(=c), through isotropic 
radiation the reciprocal of the mean free path for collisions with photons is given by 

1 f°° , n(e) /- s max(e,E 7 ) 

2~— / de — — 1 ,/„. i v - „- — » 

8E^ p 7 Je min £ 



x 77 (£ 7 ) = 2 / de—^ ds(s-mc)a(s), (1) 



where n(e) is the differential photon number density, and a(s) is the total cross section for 
a centre of momentum frame energy squared given by s = m'i f c 4l +2eE 1 (l— f3 1 cos 9) where 9 
is the angle between the directions of the energetic particle (gamma-ray) and soft photon, 
s min = (2m e c 2 ) 2 , e min = (s min -m 2 c 4 )/[2£ 7 (l+/3 7 )], and s max (e,£ 7 ) = m 2 c 4 +2e£ 7 (l+/3 7 ). 
For photon-photon pair production by gamma-rays, we take m 7 = and (3~ = 1, and the 



photon-photon pair production cross section from [[19 . 

We have calculated the mean free path for this process for interactions of gamma-rays 
with the infrared and cosmic microwave backgrounds, and this is shown in Fig. 2. Also 



shown are results |21J based on the upper and lower models for the IR background due 



to Stecker and Malkan ||, extended to shorter wavelengths. We note that above 10 TeV 
the mean free path using the new IR intensity (thick solid curve) is almost an order of 
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magnitude lower, and this has very important consequences for the HEGRA observations 
of Markarian 501. 

We show in Fig. 3 the time- averaged flux from Markarian 501 observed by HEGRA 
using the stereoscopic system of telescopes |fL8| . We note that the spectrum smoothly 
extends to at least 20 TeV. This observation was confirmed by using the stand-alone 
telescopes of HEGRA, although with less statistical significance and with rather different 
systematics [^TJ . It is important to note that the Crab Nebula, as a steady galactic source, 
has approximately the same emission level as Markarian 501 in the high state, and has 
been continuously monitored by the HEGRA telescopes pi , 22]. The Crab spectrum is 



also shown in Fig. 3 for comparison (shifted down by a factor of 100 for clarity). The 
Crab Nebula spectrum determined extends again up to 20 TeV as a straight power-law 
with a very similar spectral index, and this suggests reliable detection of photons up to 
20 TeV. For the following discussion, an understanding of the absolute energy calibration 
is of considerable importance. This was investigated in great detail by Monte Carlo 
methods, and in addition, checked by absolute measurement of the cosmic ray proton flux 
using the HEGRA telescope system in about the same energy range as in the Crab and 
Markarian 501 observations. In this energy range the cosmic ray proton flux is rather well 
measured by direct methods, and a compilation of data can be found in ref. [^]. The 
absolute energy scale of HEGRA was then estimated to have an uncertainty of 15% [24]. 
In addition, there is data available on Markarian 501 from the Whipple p5| and CAT |26| 
telescopes, and this is in very good agreement with the HEGRA data but does not reach 
as high in energy, mainly because of their lower exposure. 

We have applied the mean free path in the infrared background radiation plotted in 
Fig. 2 to "correct" the observed flux from Markarian 501 assuming a Hubble constant of 
H = 65 km s _1 Mpc -1 , and this is also shown in Fig. 3. We note that this implies a 
dramatic (by several orders of magnitude) increase in luminosity at the source (see right 
hand scale). Even taking into account relativistic beaming, with a typical Doppler factor 
of 20, this would imply an extraordinarily high luminosity for an active galaxy of the BL 
Lac class of which Markarian 501 is a member. The large correction factor is a direct 
consequence of the new IR data at 60 and 100 microns and the implied shift of the "dust 
peak" to somewhat shorter wavelengths. It is worthwhile pointing out that the correction 
factor depends on the value of the Hubble constant. If it were larger than that we have 
assumed, then the distance to the source would be smaller, and the absorption less, but 
unacceptably high values of H would be required to give a reasonable source luminosity. 
In fact, the problem may be even more severe if future infrared measurements would fill 
in the "valley" between 5 and 40 microns giving an IR background such as shown by the 
short dashed lines in Figure 1. In this case, the mean free path would be as shown by 
the short dashed curves in Figure 2, with the result that correction for absorption would 
cause an even more dramatic rise of the source spectrum above a few TeV. 
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3 Immediate consequences 



The source spectrum of Markarian 501, i.e. the observed spectrum corrected for prop- 
agation from the source, is in gross disagreement with expectations of the Synchrotron 
self-Compton (SSC) model routinely used to interpret the observed spectra of active galax- 
ies. In this model, accelerated electrons produce most of the observed radio through X-ray 
emission as synchrotron radiation, and the same population Compton scatters synchrotron 
photons to gamma-ray energies. The apparent success of the SSC model depends critically 
on which intergalactic absorption is used, as the observed sources are at the gamma-ray 
horizon 
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Now that the most recent determinations of the diffuse IR background are 
high, this makes it impossible to fit the gamma-ray observations using the SSC model. 
For example, Guy et al. took a relatively high level for the IR background, but still 
a factor of 2 short of recent determination. Petry et al. |29] and Bednarek and Protheroe 
20H used the IR background at a level of that predicted by Malkan & Stecker H, while 



Sambruna et al. |30[ assume a specific shape of the IR background in order to keep the 
TeV data compatible with the SSC model. 

The corrected TeV data (Fig. 3) have a very prominent upturn above 10 TeV which 
cannot be explained by SSC models. Such an upturn may occur naturally in some proton 
blazar models , particularly those which predict a very flat or nearly monoenergetic 
spectrum of synchrotron radiating protons [[3^, [53]. However, even with these models it 
appears impossible to obtain the extreme turn-up in the TeV source spectrum inferred if 
the new IR data are used to correct for propagation. 

In principle the turn up could be a pile-up caused in a pair-Compton cascade during 
propagation through the IR background fll3l . However, for this to be the cause the 



intergalactic magnetic field would have to be extremely small otherwise the pile-up would 
be distributed as a halo [34|, and any time- variability, such as the observed flaring on 



time scales of hours ET| would be washed out. 



4 Discussion: What is the solution? 

Many possible solutions have been mentioned by Finkbeiner et al. || , and here we provide 
a quantitative discussion of some of them. The simplest solution to the problem would be 
if the recent determinations of the universal infrared background were an overestimate. 
This is not completely ruled out, and they are not incompatible with star formation and 
dust-heating in the early Universe. In a similar manner, a shift in the energy scale of the 
HEGRA data downwards by a factor of two would have a similar effect. This appears 
rather extreme, and not very likely. If both the infrared data and the HEGRA TeV data 
are confirmed, then we may be forced to consider more radical possibilities. 

An interesting idea based on known physics was proposed by Harwit, Protheroe and 



Biermann |35| . It considers a TeV Cherenkov event as being due to a coherent super- 
position of a number of lower energy photons simultaneously arriving at the top of the 
atmosphere, and masquerading as a single TeV photon with the sum of the energies. Such 
Bose-Einstein condensates would suffer losses on a distance scale comparable to that of 
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the mean free path in the infrared radiation of the individual photons of the condensate. 
This idea can be tested on the basis of extensive simulations, and by detailed comparison 
of Cherenkov images of gamma ray showers from the Crab nebula and Markarian 501. 
Such analysis is underway. [] Even for a Bose-Einstein condensate with, on average, as 
few as two arriving photons, the effect on the intergalactic absorption of 20 TeV pho- 
tons from Markarian 501 would be dramatic. We show in Fig. 4 the result of correcting 
for intergalactic absorption for this case in an approximate way (plotted as asterisks). 
In calculating this, we assumed for simplicity that each event recorded by the telescope 
consisted of a Bose-Einstein condensate of two photons, and we then calculated the av- 
erage number Nq of photons of energy E 1 /2 that would have been emitted for each pair 
of photons arriving with total energy E 7 . The energy of the emitted bunch was then 
on average NqE^/2, and the correction factor used was therefore N /2. Because of the 
change in occupation number, and hence energy, of the bunch of photons during propaga- 
tion, the precise meaning of the corrected spectrum (E 2 F(E y )) is ambiguous, and caution 
should be applied when interpreting it. Nevertheless, integrating the corrected luminosity 
(E 2 L(E y )) over ln(i? 7 ) will give the correct total emitted power. Thus, we note that this 
process is a viable mechanism for solving the problem, and has not yet been ruled out by 
the observations. 

Another possibility which has been suggested is violation of Lorentz invariance 0, [37], 
38, |39|. Solutions of that kind have also been invoked [0] to explain the observed high en- 



ergy particles above the GZK J|l], Q cut-off. Amelino-Camelia et al. |37] suggested that 
the velocity of propagation of energetic particles is modified in vacuo due to microscopic 
quantum fluctuations occurring on scales of the order of the Planck length. Their disper- 
sion relation was used by Kifune [38 to show that intergalactic space will be much more 
transparent to TeV photons. We follow his treatment, in which energetic photons have 
velocity /3 7 c = (1 — £ y Ey/E )c, and momentum p given by p 2 c 2 = (E 2 + ^ 7 E^/E ) from 
which we can define an effective mass m 7 for high energy photons by m 2 c A = —^E^/Eq. 
Taking £ 7 = 1 and E = 1.2 x 10 19 GeV (Planck mass) |36|, |7|> and using /3 7 and m 2 
in Equation 1, we get the result shown by the thick chain line in Figure 2. One should 
be cautious, however, as Equation 1 was derived assuming s to be Lorentz invariant. 
Nevertheless, our result gives an indication of the likely effect of Lorentz Invariance vio- 
lation. Applying our result to correct the HEGRA 1997 Markarian 501 data we obtain 
the spectrum shown by the diamonds in Fig. 4, and find that indeed the absorption in 
intergalactic space is sufficiently reduced to solve the problem, resulting in a spectrum 
that is largely consistent with both leptonic and hadronic models. 

However, it is remarkable that an acceptable source spectrum for Markarian 501 is 
obtained simply by choosing the Planck mass as the energy scale for violation of Lorentz 
Invariance. Limits for this energy scale have been obtained previously by looking for 



time delays in variable sources such as Gamma Ray Bursts [f£|] and Markarian 421 [ fHfl . 
By observing instead the energy spectrum of distant sources like Markarian 501 above 
10 TeV, we may also have a sensitive tool for exploring the energy scale of Lorentz 



1 After submission of this paper a preprint by the HEGRA Collaboration appeared ( astro- ph / 000609^ ) 
which renders the possibility of a Bose-Einstein condensate explanation as being rather unlikely. 
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Invariance violation. 

New data on Markarian 501 and Markarian 421 relevant to the issue discussed in this 
paper are soon to come as both sources are observed to be in a high flaring state in the 
ongoing observation period (May 2000). Detecting TeV gamma-rays from a more distant 
BL Lac object (e.g. z =0.03-0.09) would provide for a much needed additional constraint 
on the phenomenon observed at Markarian 501. 
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Figure 1: Infrared background radiation field. Recent determinations from DIRBE || 
(diamonds), ISOCAM[0] (square), DIRBE || (triangles), and FIRAS[| (crosses) are com- 
pared with the models of Malkan and Steckerg (dot-dash curves; the dot-dot-dot dash 
curves are an extrapolation of the Malkan and Stecker models used in ref. p0 |. The lower 
limits are from ref. |TD[. In this paper we model the IR background by the solid curve; 



thick dashed curve is allowed if the ISOCAM point is considered a lower limit. 
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Figure 2: Mean free path for photon-photon pair production in the infrared-microwave 
background radiation. The curves correspond to those in Fig. 1 except that the effect of 
Lorentz Invariance violation discussed in Section 4 is shown by the long dashed curve. 



10 



CM 



> 

CD 

E- 



HO 
O 



i r 

o r 
-l [ 

-3 : : 

-4 r 

-? i- 

-8 : : 
-9 i- 
- 10 L 



-l 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

Mkn 501, 2=0.0336 
HEGRA 1997 average 



corrected 



++++++++++ 



1 
log(VTeV) 



-. 50 „ 




48 *° 



CD 



44 kT 

cm 
o 



-.42 



Figure 3: The time-averaged spectrum of gamma-rays from Markarian 501 observed in 
1997 JT8[ is compared with the spectrum of the Crab Nebula observed in 1997 8 pT], 



The spectrum of Markarian 501 after correction for absorption in the infrared background 
is also shown assuming H = 65 km s" 1 Mpc -1 . The right hand scale shows the luminosity 
for Markarian 501. 
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Figure 4: The time-averaged spectrum of gamma-rays from Markarian 501 observed by 
HEGRA in 1997 is plotted and compared with the spectrum after correction assuming that 
each observed "gamma ray" is a coherent superposition of two photons (Superposition), 
or that Lorentz Invariance is violated (LI violation); upper limits have been omitted for 
clarity from the "corrected" data. 



12 



